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ABSTRACT 

We investigate the FIR properties of a sample of BCDs observed by AKARI. By utilizing the 
data at wavelengths of A = 65 /im, 90 /xm, and 140 /im, we find that the FIR colours of the 
BCDs are located at the natural high-temperature extension of those of the Milky Way and 
the Magellanic Clouds. This implies that the optical properties of dust in BCDs are similar 
to those in the MiUcy Way. Indeed, we explain the FIR colours by assuming the same grain 
optical properties, which may be appropriate for amorphous dust grains, and the same size 
distribution as those adopted for the Milky Way dust. Since both interstellar radiation field 
and dust optical depth affect the dust temperature, it is difficult to distinguish which of these 
two physical properties is responsible for the change of FIR colours. Then, in order to exam- 
ine if the dust optical depth plays an important role in determining the dust temperature, we 
investigate the correlation between FIR colour (dust temperature) and dust-to-gas ratio. We 
find that the dust temperature tends to be high as the dust-to-gas ratio decreases but that this 
trend cannot be explained by the effect of dust optical depth. Rather, it indicates a correlation 
between dust-to-gas ratio and interstellar radiation field. Although the metallicity may also 
play a role in this correlation, we suggest that the dust optical depth could regulate the star 
formation activities, which govern the interstellar radiation field. We also mention the im- 
portance of submillimetre data in tracing the emission from highly shielded low-temperature 
dust. 
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1 INTRODUCTION 

The far-infrared (FIR) emissio n from galaxies i s often used to trace 
the star formation activities jKennicuttl [1993: llnoue et al.l [2OOO'; 
llgles ias-Paramo et al. 2003). Dust grains are the source of FIR 
emission, and the strong connection between FIR luminosity and 
star formation rate can be explained if the ultraviolet (UV) light 
from massive stars is the dominant source of dust heating. The FIR 
spectral energy distribution (SED) of dust grains reflects various 
information on the grains themselves and on the sources of grain 
heating (e.g.lTakaei. Vansevicius. & Arimotcll2003l : iTakeuchi et al.l 



I2OO5I : iDopita et alj I2OO5I) . Dust temperature is determined by the 
intensity of interstellar radiation field (ISRF) and the optical prop- 
erties of dust (i.e. how it absorbs and emits light). Since dust grains 
absorb UV light efficiently, the FIR luminosi ty and the dust te m- 
perature mostly reflect the UV radiation field jSuat & Xulll996t) . 

Observationally, dust temperature can be estimated from FIR 
colour, which is defined as the flux ratio between two FIR wave- 
lengths. If more than two (^ 3) bands are available in FIR, we 
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can take two independent FIR colours and examine a FIR colour- 
colour relation. A colour-colour relation, because of increased in- 
formation compared with a single FIR colour, enables us to obtain 
not only the dust t emperature but also the wavelength dependence 
of dust emissivity. iNagata et al] l l2002h , adopting 60 /xm, 100 /im, 
and 140 /xm as three wavelengths, show that there is a tight re- 
lation in the FIR colour-colour diagram. This tightness implies a 
common wavelength dependence of FIR e missivity among v arious 
galaxies. Their work has been developed bv lHibi et alj ( |2006|) . who 
show that a tight FIR colour-colour relation found for the Milky 
Way dust emission, called "main correlation", is also consistent 
with the FIR colour-colour relation of a sample of nearby galax- 
ies (see also ISiJcon et alj |2004 lOnaka et ai]|2007l) . As stated by 
iHibi et akl ( l2006h . this implies that the optical properties of grains 
in FIR are common among the Milky Way and nearby galaxies. 

Since the FIR colour-colour relation can be used to con- 
strain the optical properties of dust grains, the FIR colour- 
colour relation of galaxies with various evolutionary stages pro- 
vides pieces of information on the grain evolution in galactic 
environments. In particular, the dominant production source of 
dust grain could be different in early phase of galaxy evolution 
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jTodini & Ferrarj|200ll:lNo zawa et alj l2003l : iMaiolino et ai]|2004 
TSchneide r. Ferrara, & Salv aterra Moreover, the effects of in- 

terstellar processing of dust grains, especially accretion of heavy 
elements onto dust grains, coagulation, and shattering, depend on 
metallicity (or dust abundance). Thus, it is probable that the grain 
properties evolve as galaxies evolve. 

Although it is difficult to obtain the rest-frame FIR data of dis- 
tant galaxies in an early phase of galaxy evolution, there are nearby 
possible "templates" of primeval galaxies, blue compact dwarf 
galaxies (BCDs). Indeed, BCDs have o n-going star formation 
in metal-poor and gas-rich enviro nments jSargent&Seari3ll970l : 
Ivan Zee. Skillman. & SalzeJl998h . which could be similar to those 
in high-redshift star-forming ga laxies. Moreover, BC Ds harbour an 
appreciable amount of dust (e.g. lThuan. Sauvage, & Madden. 19990 
and can be used t o investigate the dust properties in chemically un- 
evolved galaxies jTakeuchi et Zll2005h . 

It is possible to investigate the metallicity dependence of dust 
prope rties by sampling BCDs with a variety of metallicity. Re- 
cently lEngeSrachTet^] I2OO8) have examined the metallicity de- 
pendence of dust emission in various wavelengths by using Spitzer 
data. In FIR, they have used Multiband Imaging Photometer for 
Spitzer (MIPS) 70 pm and 160 bands, and have shown that 
there is a correlation between dust temperature derived from these 
two bands and metallicity. This indicates the importance of studies 
on dust emission in a wid e metallicity rang e. 

The AKARI satellite (Muraka mi et alj|2007 ') provides us with 
a good opportunity to study FIR colour-colour relations, since Far- 
Infrared Surveyor (FIS) on A KARI has four bands in FIR (65 fim, 
90 Atm , 140 fim, and 160 fim) ( Murakami e t al.'2Q07':'Ka wada et"al] 
l2007h . Indeed, seven of the eight BCDs in Hirashita et alj j2008l 
hereafter H08) are detected at 65 ^m, 90 fim, and 140 pm. This in- 
dicates that it is possible to study the FIR colour-colour relation of 
BCDs by using AKARI data. In this paper, we add four more BCDs 
available in the AKARI archive and examine the FIR colour-colour 
relation of BCDs. Then, we extract information on what determines 
or regulates the FIR emission in metal-poor environments with the 
aid of theoretical models for dust emission. 

This paper is organized as follows. First, in Section[2l we de- 
scribe the models of dust emission with a simple radiative trans- 
fer recipe. Then, in Section [3] we explain the data analysis of the 
BCD sample observed by AKARI. In Section |4] we overview the 
results of the model calculations in comparison with the observa- 
tional data. In Section [5] we discuss our results, focusing on the 
relation between FIR colour and dust content. Finally, the conclu- 
sion is presented in Section|6] 



2 FIR SED MODEL 

We adopt the theoretical framework of iDraine & Lil i200lh to cal- 
culate the SED of dust emission in FIR. The framewo rk has al- 
ready been described in lHirashita. Hibi. & Shibail ( l2007l hereafter 
HHS07), but we modify it to treat BCDs in this paper. Moreover, as 
described in the last paragraph in Sect. 12.11 we newly include the 
effect of radiative transfer in a similar way to Galliano et al. ( 2003). 
We overview our framework, focusing on the change from HHS07. 

2.1 Interstellar radiation field 

The stell ar SEDs of BCDs are generally har d er than those of spira l 
galaxies dBoselU. Gavazzi. & Sanvitoll2003h . iMadden et ai] ( l2006h 
also show a hard stellar radiation field from the observation of 



mid-infrared emission lines in low-metallicity star-forming galax- 
ies. Thus, we change the SED of ISRF, and adopt the following 
fitting formula according to the averaged SED of the BCD sample 
in ,Boselli et a l. ( 2003): 

(A^,„ ^ 0.0912), 
3.02 X W-*\-^'^° 

(0.0912 < < 0.25), 



inXJ 



(0) 



(1) 



1.55 X lO"'' X 10"'=' A;^"'''^ 

(0.25 < X^m < 0.37), 
2.39X1^^ (0.37 < A^„, < 0.4), 
inX[WiBx{n) +W2Bx{T2)] 

(0.4 < A,„n), 

where J^^ is the ISRF intensity whose normalization is determined 
so that the inte nsity at wavelengt h A = 0.2 /im (almost the centre of 
the UV range; lBuat&Xulll99d) is equal to the solar neighborhood 
ISRF estimated bv lMathis, Mezger, & Panagial ( Il983l) (47rA J^"' is 
expressed in units of erg cm~'^ s~^), A^m is the wavelength in units 
of ^im,BA(T) is the Planck function, (Ti, T2) = (1000, 4000 K), 
and {Wi, W2) = (1.61 X 10"l^ 2.90 x 10"'''). 

In this paper, we newly include the effect of dust extinction on 
the ISRF. We simply assume that the ISRF intensity, J\ , at optical 
depth Tx is 



(2) 



where x is the scaling factor of the ISRF relative to J^*^' and is 
assumed to be independent of wavelength. For the wavelength de- 
pend ence of ta, we assume the M i lky W ay extinction curve taken 
from lCardeUi. Clavton. & MathisI ( Il989h with Rv = 3.1. Thus, 
the extinction at any wavelength A, A\ (in units of magnitude), can 
be determined if we set a value of Av (extinction at V band). The 
optical depth can be related with the extinction as ta = Av/ 1-086. 
The following results does not change drastically if we assume 
an extinction curve appropriate for the Large Magellanic Clouds 
(LMC) or the Small Magellanic Cloud (SMC), since the detailed 
shape of UV-optical SED is not important in this paper Strictly 
speaking, we treat scattering as effective absorption, and this treat- 
ment would overestimate the absorbed energy. Thus, we interpret 
Av given here as effective absorption optical depth, so that the total 
energy absorbed by dust in our formulation is equal to the energy 
that would be absorbed if scattering were properly treated. 

iGalliano et al.l ( l2003h also adopted a similar treatment, which 
is valid if dust is distributed in a thin shell surrounding the stars. 
Such a "screen" geometry enhances the effect of dust extinction 
(i.e., radiative transfer) in comparison with a mixed geometry be- 
tween stars and dust. Thus, our models are suitable to examine the 
extent to which the radiative transfer effects could affect the FIR 
SEDs. 



2.2 Grain properties 

We consider silicate and graphite grains in this paper Since we 
are interested in the wavelength range appropriate for the AKARI 
FIS bands (A ~ 50-180 pm; Section l3Tt . we neglect polycyclic 
aromatic hydrocarbons (PAHs) , which contribute significantly to 
the emission at A < 20 / xm (Deser t. Boulanger. & Pugetlll990l : 
iDwek et al.lll997l : lDraine & Li 2001). 

We assume a grain to be spherical with a radius of a. The 
absorption cross section of the grain is expressed as 7ra^(5abs(A), 
where Qabs(A) is called absorption efficiency. We adopt the op- 
tical constants of astronomical silicate and graphite for A < 
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lOOum from lOraine & Le3 ( 1 19841) and IWeingartner & Drain^ 
( 12001) and calculate the absor ption efficiency by using Mie the- 
ory dBohren &Huffmanl 11983'). For the absorption efficiency at 
A > 1 00 iim, we adopt a functional form proposed by , Reach et al.i 
( Il995h : 



Qabs(A) 



(A/Ao) 



[l + (Ai/A)«]i/6 



(3) 



which behaves like P — 1 ((3abs(A) oc A"'', where /3 is called 
emissivity index) for A <C Ai and /9 = 2 at A 2> Ai (see also 
iBianchi. D avies. & Alton" 1999). We assume that Ai = 200 /im, 
following Reach et al. ( 1995). We set the value of Ao so that the 
continuity of Qabs at 100 is satisfied. As shown in HHS07, 
the wavelength dependence of Qabs assumed here explains the FIR 
colour-colour relations of the Milky Way, the LMC and the SMC. 
Although the absorption efficiency at A > 100 /im is modified, 
we call the grain species "silicate" and "graphite" according to the 
adopted absorption efficiency at A < 100 ^m. 

The number density of grains with sizes between a and a + 
da is denoted as Ui (a) da, where the subscript i denotes a grain 
species (silicate or graphite). Since there is little knowledge on the 
grain size distribution in BCDs, we simply assume a power-law 
form for ni(a): 



ni{a) = da ^ (amin < a ^ amax) , 



(4) 



where amin and amax are the upper and lower cutoffs of grain 
size, respectively, and d is the normalizing constant. We assume 
K = 3.5, a min = 3.5 A and amax = 0-25 am for both graphit e 
and silicate ( Mathi s. Rumpl. Nordsiec3ll977l : ILi & Drainell200ll) . 
As long as we treat the FIR colours of a single species (Section 
12.3b , the constant d cancels out. Because the dust composition in 
BCDs is still unclear, we consider silicate and graphite separately 
as two possible dust species. 

Following HHS07, we adopt the multi-dimensional Debye 
models for the heat capacities of grains ( Draine & Li 2001i) . The 
physical parameters necessary to characterize the heat capacities 
can be found in HHS07. 



2.3 Calculation of FIR colours 

The FIR intensity (per frequency per hydrogen nucleus per 

solid angle) of the emission from dust species i at a wavelength A 
can be estimated as 



da — ni(a)Tva Qabs(A) 



dP 

dTB^iT)^, 



(5) 



where Bv{T) is the Planck function at frequency v and temper- 
ature r, and APi/AT is the temperature distribution function of 
the grains under the ISRF evaluated in Section lTTI The calculation 
method of APi/AT is based on the concept of stoc hastic heating, 
and we follow the formalism of lOraine & Lil j200lh (see their sec- 
tion 4). We assume that the FIR radiation is optically thin. Thus, 
after the radiative transfer effect, the intensity at A per frequency 
per hydrogen nucleus per solid angle, /^(A), becomes 



X As , 



(6) 



where s is the path length, which is proportional to Av- Here the 
colour at wavelengths of Ai and A2 is defined as I^i^i) / Iv{^2), 
which is called Ai — A2 colour and denoted as (Ai / A2)ci (we omit 



^m in this expression; e.g. (60/100)ci and (140/100)ci). Since we 
only consider colours (i.e. flux ratios), the proportionality constant 
between s and Ay does not matter in this paper. Thus, practically, 
we can replace s with Av- 

Since we calculate colours in this paper, we can de- 
fine the colour with Av = as the optically thin extreme 
(limAv,-,o(Ai/A2)ci). In other words, the FIR colour calculated 
with Av ~ (in fact Av 0) represents the case where the dust 
is illuminated by a radiation field without extinction. In the follow- 
ing, Av ~ means that the dust is illuminated by a radiation field 
without extinction, not that there is no dust. 



3 DATA 

3.1 Analysis of the AKARI data 

H08 analysed eight BCDs observed by FIS onboard AKARI with 
four photometric bands of N60, WIDE-S, WIDE-L, and N160, 
whose central wavelengths are 65, 90, 140, and 160 nm with effec- 
tive band widths of AA = 21.7, 37.9, 52.4, and 35.1 fim, respec- 
tively dKawada et alj|20()7l) . The measured FWHMs of the point 
spread function are 37 ± 1", 39 ± 1", 58 ± 3", and 61 ± 4", re- 
spectively, for the above four bands. We exclude Mrk 36 from the 
sample, since it is not clearly detected at WIDE-L (140 /im). The 
other seven sample BCDs in H08 are listed in the lower part of 
TableE 

In this paper, we add four BCDs available in the AKARI 
archivtQ: UM 420, Mrk 59, Mrk 487, and SBS 1319-1-579 as shown 
in the upper part of Table [T] We selected these galaxies among the 
BCDs observed by FIS after investigating if they were detected at 
N60 (65 /im), WIDE-S (90 ^im), and WIDE-L (140 fim). The ob- 
serving mode (FISOl), the scan speed (8" s~^), and the reset inter- 
val (2 s) are the same as those adopted for the sample in H08. The 
method of data reduction and analysis is the same as that in H08, 
and the measured fluxes are listed in Table [T] Below we briefly re- 
view the data reduction and analysis (see H08 for details). 

The raw data were reduced by using the FIS Slow Scan Tool 
(version 20070914)0 Because the detector response is largely af- 
fected by a hit of high energy ionizing particle, we used the lo- 
cal flat; that is, we corrected the detector sensitivities by assuming 
uniformity of the sky brightness jVerdugo, Yamamura. & PearsonI 
l200l) FI The errors caused by the smoothing procedures are well 
within 10% for ^60 (65 ^im) and WIDE-S (90 ^im), 20% for WIDE- 
L (140 /im), and 25% for N160 (160 /im). These error levels are 
consistent to or somewhat larger than the background fluctuation 
mainly caused by a hit of high energy ionizing particle, which 
implies that the background fluctuation is a major component in 
the errors. The above values are adopted for the errors, but we 
adopt larger errors for some galaxies which have larger background 
noises. At ^60 (160 /im), some objects are not detected. For them, 
we adopt 3 times the background uncertainty as an upper limit. 
Since the sample BCDs are compact enough to be treated as point 



^ http ://darts . isas .j axa.jp/a.stro/akari/. 

^ http://www.ir.isas.ac.jp/ASTRO-F/Observation/. 

We refeiTed to AKARI FIS Data User Manual Version 2 for the data 



reduction and analysis jhttp://www.ir.isas. ac.jp/AST RO-F/Observation/ 1. 
When we ran the pipehne command ss_run_ss, we applied options 
/local, /smooth, width.f ilter=80: The flat field was built from 
the observed sky, and boxcai' smoothing with a filter width of 80 s in the 
time series data was applied to remove remaining background offsets among 
the pixels. 
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Table 1. Measured fluxes ("<" indicates the 3 cr upper limit). Tlie upper four objects are newly 
analysed in this paper, while the lower seven are the sample in Hirashita et al. ( 2008). 



Name 


N60 


WIDE-S 


WIDE-L 


N160 




(A = 65 fim) 


(A = 90 fim) 


(A = 140 fim) 


(A = 160 nm) 


UM420 


0.46 ± 0.05 Jy 


0.37 ±0.04 Jy 


0.14 ±0.03 Jy 


< 0.36 Jy 


Mrk59 


2.3±0.2Jy 


2.1 ±0.2 Jy 


1.4±0.3Jy 


1.2±0.3Jy 


Mrk 487 


0.45 ± 0.05 Jy 


0.43 ± 0.04 Jy 


0.21 ±0.06 Jy 


< 0.69 Jy 


SBS 1319+579 


0.38 ± 0.04 Jy 


0.29 ± 0.03 Jy 


0.28 ±0.06 Jy 


< 0.59 Jy 


11 Zw 40 


6.9±0.7Jy 


6.6±0.7Jy 


3.7±0.7Jy 


3.4±0.9Jy 


Mrk 7 


0.83 ± 0.23 Jy 


0.84 ± 0.08 Jy 


1.1 ± 0.2 Jy 


1.0±0.3Jy 


Mrk 71 


3.7±0.4Jy 


3.1 ±0.3 Jy 


2.2 ± 0.4 Jy 


< 3.0 Jy 


UM439 


0.42 ±0.04 Jy 


0.40 ±0.04 Jy 


0.31 ± 0.06 Jy 


< 0.45 Jy 


UM533 


0.54 ±0.05 Jy 


0.61 ±0.06 Jy 


0.54 ± 0.11 Jy 


< 0.68 Jy 


II Zw 70 


0.89 ±0.09 Jy 


0.80 ±0.24 Jy 


0.75 ± 0.15 Jy 


0.62 ± 0.16 Jy 


II Zw 71 


0.32 ±0.03 Jy 


0.54 ±0.05 Jy 


1.3 ±0.26 Jy 


1.3±0.3Jy 



sources with the resoluti on of FIS, w e follow the point-source pho- 
tometry described in Ver dugo et al 

Finally the same colour correction as in H08 was adopted for 
the WIDE-L (140 /xm) fluxes: the colour correction factor was as- 
sumed to be 0.93 (the flux was divided by this factor). We did not 
apply colour correction to the N60 (65 fim), N160 (160 /im), and 
WIDE-S (90 fim) fluxes, following H08. For these three bands, the 
uncertainty caused by not applying colour correction is smaller than 
the errors put in Table [T] 



3.2 FIR colours 

The FIR colours are defined by the flux ratio between two bands 
(Section [23J. HHS07 investigated the FIR colours of the Milky 
Way and the Magellanic Clouds by using the Zodi-Subtracted Mis- 
sion Average (ZSMA) taken by the Diffuse Infrared Background 
Experiment (DIRBE) of the Cosmic Background Explorer (CO BE). 
Detail s of the observational data analysis can be found in Hibi et al.l 
( l2006h . who adopted DIRBE bands of 60 ^m, 100 fim, and 140 fim 
and two colours, 60 /im-100 ^m colour, (60/100)ci, and 140 fLm- 
100 fim colour, (140/100)ci. 

We can also take two colours for our BCD sample. However, 
the wavelengths of the FIS bands are slightly different from the 
DIRBE bands. Thus, we apply the following corrections. First, we 
fit Av^B,^ (Td) {A is a constant and (3 is the emissivity index) to the 
65 fim and 90 fim fluxes derive the dust temperature Td. We denote 
the dust temperature obtained in this way as Td((65/90)ci, /3). We 
adopt 13 = 1 and 2. Then, we estimate the 60 fim flux by eval- 
uating Aul^B^iTa) at A = 60 fim under the values of A and 
Td((65/90)ci, P) obtained above. The same procedure is applied 
for the 140 ^m and 90 fim fluxes to obtain the 100 fim flux. In this 
way, we can obtain (60/100)ci and (140/100)ci for the BCD sam- 
ple. In Fig.[T] we show the colour-colour relatio n for the sample. 

For comparison, we also show the data of iHibi et al.l ^20061) 
for the Milky Way (the Galactic plane with Galactic latitudes of 
|6| < 5°) and the Magellanic Clouds in Fig.fT] Hibi et al. (2006) 
found that more than 90% of the data lie on a strong correlation 
called main correlation, which can be fitted as 

(140/100)ci = 0.65(60/100);;;°-^* . (7) 

The main correlation also explains the FIR c olours of Galactic high 
latitu des (|&| > 5°), the LMC, and the SMC faibi et al.ll2006l : lHibil 
l2006l : HHS07). In the Galactic plane, there is another correlation 




0.5 1 5 

(140/100),, 



Figure 1. Far-infrared colour-colour diagram showing the relation between 
60 fim-lOQ fan colour and 140 /im-100 fim colour. The filled squares with 
error bars indicate the current BCD sample, while the crosses and open dia- 
mon ds show the DIRB E data of the LMC and the SMC, respectively, taken 
from Hibi et al ] <2006l) . We also presen t the Milky Way data (with Galac- 
tic latitudes of < 5°) adopted from lHibi et al 1 i2006l) by the contours, 
which show the levels where 50, 80, 90, 95, and 99 per cent of the data are 
contained. The dashed line and the dotted lines present the fittings for the 
main coiTelation (equation |7) and the subcorrelation (equation [s), respec- 
tively. 

sequence, called subcorrelation in lHibietalJ ( l200^ : 

(140/100)ci = 0.93(60/100)~i° '''^ . (8) 

The subcorrelation is not seen in the high Galactic latitud es 
I2OO6I : HHS07). This supports the idea of Hibi et aP j2006h that the 
subcorrelation is produced by a contamination of high-ISRF re- 
gions, which tend to reside in the Galactic plane. 

It is interesting that not only the LMC and the SMC but also 
the current BCD sample has consistent FIR colours to the main cor- 
relation or the subcorrelation (Fig.[T]l. This implies that the wave- 
length dependence of the FIR emissivity is not different among the 
BCDs, the LMC, the SMC, and the Milky Way. In the following 
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Figure 2. Far-infrai'ed colour-colour relations for silicate and graphite (as- 
terisks and open diamonds, respectively). The different points correspond 
to X = 1, 3, 10, 30, 100, and 300 from high to low [140/100]. We adopt 
Av = (optically thin extreme) for the ISRF. The same data as shown 
in Fig. [T] as well the the main conelation and subcorrelation (dashed and 
dotted lines, respectively) are also presented. 



section, we examine if the FIR colours of BCDs can really be re- 
produced with the emission properties adopted by HHS07, who ex- 
plained the FIR colours of the Milky Way, the LMC, and the SMC 
(Section|2ll. 



4 RESULTS 

4.1 Interstellar radiation field 

We present the dependence of FIR colours on ISRF. Here we adopt 
Av = (i.e. the optically thin extreme for the ISRF; Section [23t 
to concentrate only on the effects of ISRF intensity. In Fig. [21 we 
show the FIR colour-colour relation ((60/100)ci and (140/100)ci) 
calculated for silicate and graphite. We observe that the FIR colours 
of the BCD sample can be reproduced with x ~ 30-300 except 
for the lowest data point (II Zw 71). This is consistent with the 
conclusion in H08 that the ISRF (especially the UV ISRF, which 
contributes most to the dust heating) in BCDs can be 100 times 
higher than the Galactic ISRF in the solar neighbourhood. 

Our results not only confirm the warm IRAS {Infrar ed Astro- 
nomic al Satellite) 60 /im-100 colours of BCDs (.Hoffman et al.l 
Il989t) . but also show high temperatures of dust contributing to the 
emission at A > 100 /xm. A few Virgo Cluster BCDs observed by 
ISO {Infrared S pace Observatory) sh ow relatively cold 170 /im- 
100 jim colours i lPopescu et a"i]|2002l) . indicating that the presence 
of warm dust is not necessarily the case at least in cluster environ- 
ments. 

We also find that the calculated FIR colours cover the area 
between the main correlation and the subcorrelation. With x < 
10, the colours trace the main correlation, which is appropriate for 
the Milky Way, the LMC and the SMC. On the other hand, the 
FIR c olours are rather similar to the subcorrelation for x ^ 100. 
While lHibi et"ZI ( |2006|) and HHS07 claimed that the subcorrelation 



can be reproduced with dust emission with multiple temperature 
components, we see here that the upper part of the subcorrelation 
can be explained with a single temperature. 

4.2 Radiative transfer effects 

If radiative transfer effects are considered, FIR emission comes 
from dust at different optical depths, where the ISRF intensities are 
different because of dust extinction. As modeled in Section [23l the 
FIR int ensity is the sum of all contributions from different optical 
depths. iHibi et alj {iQQ&i expect that radiation transfer effects shift 
the FIR colours toward the subcorrelation on the colour-colour di- 
agram. 

In Fig. [3] we show the FIR colours for x = 3, 30, and 300 
with various values of Av {Av = 0, 0.1, 0.2, 0.5, 1, and 20. The 
FIR colours change little for Av > 2, since the emission is negli- 
gible from such large optical depths in the wavelength r ange of in- 
terest because of low dust temperatures. As expected bv lHibi et al.l 
( l2006l) . there is a slight trend that the colours move toward the sub- 
correlation for X J; 30. However, we observe that the radiative 
transfer effects shift the FIR colours along the main correlation or 
the subcorrelation for x ^ 30. This indicates that the radiative 
transfer effects are hard to be distinguished from the change of x 
on the colour-colour diagram. 

How can we distinguish the variation of Ay from that of x? 
We expect that if the change of Ay controls the colour variation 
there may be a correlation between dust content and FIR colours. 
Below we further discuss the observed FIR colours in BCDs in 
terms of the dust content. 



5 DISCUSSION 

5.1 Dust content and FIR colours 

The dust mass of the current BCD sample is estimated by the 
WIDE-S (90 /xm) and WIDE-L (140 fj.m) fluxes, which are suit- 
able to trace the total mass of large grains which have the dominant 
contribution to the total dust content ( Galliano et al. 2005)0 The 
dust mass Ma is related to the flux as (e.g. lHildebrandil 19831 : H08) 



i^vBv{Td) 



(9) 



where is the mass absorption coefficient of dust grains. The 
mass absorption coefficient is related with the absorption efficiency 
Q„ as (e.g. H08) 



4as 



(10) 



where s is the grain material density. If the grain rad ius is much 
smaller than the wavelength, Q^/a is independent of a dHildebrmdl 
Il983h . The values adopted in Section [Z2l indicate Q^/a = 1.14 x 
lO"'^ /im~^ for silicate and 1.23 x 10^^ ^mT^ for graphite at 
A = 125 ^m. By adopting typical densities for silicate and graphite 
as s = 3.3 g cm"'' and 2.25 g cm"'', respectively, we obtain the 



is considered as the 



* As stressed in Section 1231 the case of Av 
o ptically thin ext r eme fo r the ISRF. 

° iGalliano et al ] i2005l) also introduce a very cold dust component con- 
tributing to the submillimetre flux, but its abundance is sensitive to the as- 
sumed emissivity index of large grains (Section |53). Thus, we here con- 
centrate on the dust component traced in FIR and do not consider very cold 
dust. 
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Figure 3. Far-infrared colour-colour relations for (a) silicate and (b) graphite (asterisks and open diamonds, respectively). The results for x = 3, 30, and 300 
are shown. For each value of x, Ay = 0, 0.1, 0.2, 0.5, 1, and 2 correspond to the data points from low to high (140/100)ci. Note that Ay = is the optically 
thin extreme for the ISRF. The same data as shown in Fig.[T]as well the the main correlation and subcorrelation (dashed and dotted lines, respectively) are also 
presented. 



mass absorption coefficients at A = 125 fim as K125 = 25.9 and 
41.0 cm^ g^^ for silicate and graphite, respective ly. H08 derived 
K125 = 18.8 cm^ g~^ based on lHildebrand h983l) . 

Following H08, we adopt the following expression for Ki, to 
estimate the dust mass from the AKARI data: 



Ki, — K125 



A 



125 fim 



-13 



(11) 



The absorption coefficient adopted in Section |Z2l cannot be fitted 
with a single power law, but the index lies between /3 = 1 and 2. 
Thus, if we adopt the wavelength dependence described in Section 
12.21 the estimated dust mass should be between the dust masses 
obtained with {3 = 1 and {3 = 2. We prefer to adopt this simple 
power-law expression for here, since it can be directly com- 
pared with the results in H08 and the simple expression is useful 
for observational usage. 

In Table[2l we list the dust mass estimated from equation 
We determine the dust temperature Td by fitting Au^ B,^{Ti) {A 
is a constant) to the 90 and 140 /im fluxes. Note that A can 
be eliminated if we take the flux ratio of those two bands. The es- 
timated dust temperatures are listed in Table[2l where rd(/3 = 1) 
and Td{l3 = 2) denote the dust temperatures evaluated with 13=1 
and 2, respectively. For the mass absorption coefficient, we adopt 
K125 = 25.9 cm^ g~^ (i.e. the value for silicate), but the dust 
mass is simply proportional to Kj~25 if another value of K125 is 
adopted. The distance is taken from H08 for the common sample, 
and is estimated by using the Galactocentric velocity taken from the 
NASA/IPAC Extragalactic Database (NEDjfl by assuming a Hub- 
ble constant of 75 km s^^ Mpc^^ . The dust masses estimated with 
13= 1 and 2 are denoted as Afd(/3 = 1) and A/d(/9 = 2), re- 
spectively. Although the complex dependence of quantities in equa- 
tion ^ makes it hard to obtain an analytical estimate of the uncer- 
tainty in A/d, we have confirmed that the uncertainty is around a 

^ http://nedwww.ipac.caltech.edu/. 



factor of 2 by varying the measured fluxes within the errors. From 
Table|2l we observe that A/d(/3 = 1) is systematically smaller than 
Md{f3 = 2) by a factor of ~ 2.5, mainly because of higher dust 
temperature for (3 = 1. 

As an indicator of dust abundance, we take dust-to-gas 
ratio. For the gas mass, we adopt Hi gas mass Mhi esti- 
mated from Hi 21 cm emission observations. The data of A^hi 
are already compiled by H08 for the common sample (i.e. 
the latter seven BCDs in Table |2j. Among the newly added 
sample, there is no available Hi data for UM 420 and Mrk 
59, while the data for Mrk 487 and SBS 1319-1-5 7 9 are ob- 
tained from Hopkins, Schulte-Ladbeck. & Drozdovskvl l l2002h and 



iHuchtmeier et al 



ll2007h. respectively (for the latter, we used the 



formula in iLisenfeld & Ferraral 1 19981 to convert the H I flux to 
A/hi). The dust-to-gas ratio D is defined as T> = Md/Mui- We 
adopt Md{(3 = 1) for the dust mass, but the following discussion 
does not change if we adopt Md{l3 = 2). 

Fig.|4]shows the relations between FIR colour ((60/100)ci or 
(140/100)ci) and dust-to-gas ratio. We observe correlations with 
correlation coefficients r = —0.65 for the log(60/100)ci~log I? 
relation and r = 0.89 for the log(140/100)ci-log 2? relation. 
These correlations indicate that the dust temperature tends to be 
high in dust-poor objects. We present not only the sample BCDs 
but also the data of the Milky Way, the LMC and the S MC. The 
dust-t o-gas ratio of the Milky Way is assumed to be 0.006 dSpitzed 
Il978l) . while the dust-to-gas ratios of the LMC and the SMC are 
assumed to be 1/3 and 1/5 times of the value of the Milky Way ( Peil 
J992). For the FIR colours, we adopt the peak of the contours in 
Fig.[T]for the Milky Way, and the average of the data points for the 
LMC and the SMC. Those three points also follow the trend found 
by the BCDs. 

A possible interpretation of the correlation in Fig. |4] is that 
the relatively low dust temperatures in dust-rich BCDs result from 
the shielding (extinction) of stellar radiation. In order to examine 
this possibility, we calculate the relation between FIR colours and 
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Table 2. Some derived and compiled quantities. 



Name 




Td(/3 


= l)b 


Td(/3 


= 2)'' 


Md(/3 = 


= 1)- 


Md(/3 = 2)^^ 




12 + log (O/H) <= 




[Mpc] 


[K] 


[K] 


[Af©] 


[Mq] 






UM 420 


234 


97^ 


-249 
-31 


5o: 


h20 
-9 


2.5 X 


IQS 


8.7 X 10 




7.93 ± 0.05 


Mrk59 


10.9 


44 


fl5 
-7 


32 


+7 
-3 


2.2 X 


10* 


6.2 X 10* 




7.99 ±0.01 


Mrk 487 


11.0 


62" 


f 72 
-16 


4o: 


H6 
-7 


1.7 X 


10^ 


5.0 X 10^ 


8.0 X 10^ 


8.06 ± 0.04 


SBS 1319+579 


29.0 


34 


+7 
-4 


26 


+i 
-2 


6.8 X 


10* 


1.8 X 10^ 


1.7 X IQS 


8.10 ±0.01 


II Zw 40 


9.2 


53lf 


36 


-fS 

-4 


2.8 X 


10-4 


8.3 X 10* 


2.0 X 10* 


8.15 ±0.02 


Mrk 7 


42.4 


28 


+4 
-2 


23 


+3 
-1 


1.0 X 


iqs 


2.7 X 10^ 


3.6 X 10^ 


8.54 ±0.04 


Mrk 71 


3.4 


42 


fio 

-5 


31 


-fS 
-3 


3.9 X 


10^ 


1.1 X 10* 


1.2 X 10^ 


7.83 ± 0.02 


UM439 


13.1 


39 


fio 

-5 


30 


+4 
-3 


9.8 X 


10^ 


2.7 X 10* 


1.7 X 10* 


7.98 ± 0.03 


UM533 


10.9 


36 


+7 
-5 


28 


+4 
-3 


1.6 X 


10* 


4.2 X 10* 


5.8 X 10^ 


8.10 ±0.04 


II Zw 70 


17.0 


34 


fl5 
-5 


27 


+8 
-4 


5.9 X 


ID* 


1.5 X 10^ 


3.6 X 10* 


8.11 ±0.04 


IIZw71 


18.5 


22 


+2 

-2 


19 


+ 1 
-2 


7.4 X 


10^ 


1.8 X 10^ 


7.3 X 10* 


8.24 ±0.04 



Distance estimated by using the Galact ocentric velocity take n fro NED by assuming a Hubble constant of 75 km s~* Mpc~^ except 
for Mrk 71, whose distance is taken from lTolstov et aljj 19951) . 

^ Dust temperature estimated from the 140 /zm-lOO iim colour (Section fS.U . The emissivity index adopted is also indicated (/3 = 1 
or 2). 

Dust mass estimated in Section lsTI The emissivity index adopted is also indicated (/3 = 1 or 2). 

H I gas mass estimated from H I 2 1 cm line observations. For the latter seven BC Ds, the data are compiled in |Hirashitaetai]<2008h . 

The data for Mrk 487 and SBS 1 319+579 are obtained fromlHopkins et alj |2002|) and lHuchtmeier et alj j2007l) . respectively. 

° Oxygen abundance compiled in lHopkins et alj i2002h and[Hirashita et al. 1 2002 ) (and t he references therein). When the error estimate 

is not available, we put a typical uncertainty in the measurement (0.04) tShi et al .1200^ . 




Figure 4. Relation between far-infrared colour ((60/100)ci and (100/140)^1 in Panels a and b, respectively) and dust-to-gas ratio. The filled squares with 
vertical error bars represent the current BCD sample whose H I mass is available. The typical error of dust-to-gas ratio which comes from the error of estimated 
dust mass is shown by the horizontal bar. The data points for the Milky Way, the LMC and the SMC are also presented by the open diamonds. The solid and 
dotted lines show the colour variation expected from the model described by equation )13t with the optically thin extreme ( Ay = 0) for silicate and graphite, 
respectively. The correlation coefficients for the BCD data are r = —0.65 and r = 0.89 for Panels a and b, respectively. 



Av by using the models (same as Fig. |3j- The results are shown 
in Fig.|5] where we find that both the variation of (60/100)ci and 
that of (140/100)ci with a single value of x are too small to ex- 
plain the observed large diversity in these colours in Fig. |4] The 
colours are not so sensitive to the change of Ay partly because the 
strongest contribution from Av ^ 0, where the dust temperature is 
the highest, is always present. In particular, if Av > 2, the colour 
is insensitive to Ay since the dust temperature at such deep optical 
depths is too low to contribute to the total emission. Thus, the vari- 



ation of the FIR colours cannot be reproduced only with a variation 
of dust optical depth, but rather it should reflect the correlation (i) 
between dust-to-gas ratio and dust heating itself (i.e., ISRF), or (ii) 
between dust-to-gas ratio and dust properties (i.e., absorption coef- 
ficient and/or grain size distribution). 

It is difficult to survey all the possible cases for (ii), since there 
are few, if any, constraints on the dust properties in BCDs. How- 
ever, we can discuss some cases for (ii) as follows. HHS07 show 
that (140/100)ci changes only slightly along with the change of 
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Figure 5. Far-infrared colours ((60/100)ci in Panels a and b, and (100/140)ci in Panels c and d) for a foreground screen as a function of Ay. For the grain 
species we examine silicate (Panels a and c) and graphite (Panels b and d). In each panel, different lines show different ISRF intensities of x = 1, 3, 10, 30, 
100, and 300. 



grain size distribution, since it reflects the equilibrium grain tem- 
perature. Thus, it is hard to reproduce the clear correlation between 
(140/100)ci and dust-to-gas ratio only with the change of grain 
size distribution. The change of grain material as a function of 
dust-to-gas ratio may change the grain absorption coefficient, lead- 
ing to variation of dust temperature as a function of dust-to-gas 
ratio. However, we have shown that the colour-colour relation of 
the BCDs can be explained by the emissivity that also explains the 
colour-colour relation of the Milky Way. In this context, it is not 
probable that the dust property is the main driver for the change of 
FIR colours. 

The possibility (i) is interesting to investigate, since the corre- 
lation in Fig.|4]is a natural extension of the sequence of the Milky 
Way, the LMC, and the SMC. For those "nearby" thr ee galaxies, 
it is k nown that the typical ISRFs are different (e.g. IWeltv et al] 
l2006h . and indeed the FIR colour sequence can be explained by the 



variation of the ISRF l lHibi et al.|[20o3. HHS07). Thus, in the fol- 
lowing subsection, we discuss a possibility that the ISRF changes 
as a function of dust-to-gas ratio. 

5.2 Correlation between dust-to-gas ratio and ISRF? 

Here we interpret the correlation between dust-to-gas ratio and 
FIR colours as the correlation between dust-to-gas ratio and ISRF. 
Since the main heating source of dust grains is UV pho tons com- 
ing from the young stellar population jSuat & Xulll996h . the dust- 
heating ISRF reflects the spatial concentration of recent star for- 
mation activity. Thus, the correlation between dust-to-gas ratio and 
ISRF indicates that dust-poor galaxies tend to host concentrated 
star-forming regions. 

There are some theoretical models which suggest that dust 
plays an active role in determining the properties of star forma- 
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tion. Indeed lHirashita & Huri3 l l2004h show that the shielding of UV 
heating photons plays an important role in determining the strength 
of star formation: If the optical depth of dust becomes large, the UV 
heating photons are efficiently blocked and the gas efficiently cools 
to produce a favourable condition for the star formation. Thus, the 
optical depth of dust in UV is important, and a burst of star forma- 
tion could become possible if the UV optical depth exceeds one. In 
this scenario, the condition for star formation can be written as 



KVpR ~ 1 , 



(12) 



where k is the UV absorption coefficient per dust mass and R is the 
size of the star-forming region. On the other hand, we assume that 
a constant fraction of gas mass is converted to stars: L oc A/gas, 
where Mgas (~ pR^) is the gas mass, and L is the stellar UV lumi- 
nosity. Since x oc L/R^ oc M^^/R? ~ pR, we obtain x oc TJ~^ 
by using equation M2\ . Normalizing the quantities to the Milky 
Way values (i.e. x ~ ^ ^"d V = 0.006), we obtain 



V 0.006/ 



(13) 



In Fig. [4] we calculate the relation between the FIR colours 
and D by assuming equation ( 113b . The optically thin extreme (i.e. 
Av ~ 0) for the ISRF is adopted to focus on the effect of x- The 
results indeed traces the trend of the observational data. Thus, the 
above argument on dust shielding of UV photons (equation 1 12t as 
the condition for a burst of star formation is compatible with the 
data. 

Equation l ll2t implies that the size-density relation of 
star-forming region becomes roughly a relation with constant 
column density. Indeed, the observational data presented by 
iHirashita & Hunll (12009') support consta nt column density of star- 
forming regions in BCDs (i.e. p oc fl~^). |Pak et ai] (IT998) also find 
that the column density of gas is inversely proportional to the dust- 
to-gas ratio based on an analysis of photodissociation regions in the 
Milky Way, the LMC, and the SMC, which is consistent with our 
picture. 'McKee (1989) show that an equilibrium between the grav- 
itational contraction and the energy injection from stars is achieved 
at Av ~ 4-8 based on the stability argument on photoionization- 
regulated star formation, in which shielding of photoionizing pho- 
tons permits ambipolar diffusion to proceed a nd stars to form . Al- 
though the physical processes considered in iMcKed ( Il989h are 
more detailed and sophisticated, the necessity of shielding is com- 
mon with our picture expressed in equation l ll2t . 

5.3 FIR colour and metallicity 

The observational correlation between FIR colour and dust-to- 
gas ratio shown in thi s paper has confirmed the results of 
lEngelbracht et al.l ilOOA who show the correlation between 
Spitzer MIPS 70 /im-160 pm colour and metallicity. Since there 
is a correlation between dust-to-gas ratio and metallicity as shown 
by them and also by H08, the correlation with metallicity is equiv- 
alent to the correlation with dust-to-gas ratio. For the detailed dis- 
cussions on the relati on between dust-to-gas rati o and metallicity 
in BCDs, see H08 and lLisenfeld & Ferraraifl998l) . 

Fig. [6] shows the relation between FIR colour and metallicity 
(oxygen abundance in gas phase) for the BCD sample. The data 
of oxy gen abundance are compil ed in lHirashita, Taiiri, & Kamaval 
( I2OO2I) and lHopkins et al.l ( l2002h , and summarized in Table |2] We 
observe correlations with r = —0.45 for the relation between 
log(60/100)ci and 12 + log(0/H) and r = 0.53 for the rela- 
tion between log(140/100)ci and 12 + log(0/H). We also show 



the data points for the Milky Way, the LMC, and the SMC. The 
Milky Way gas metallicity is represented by the Orion Nebula 
(12 + log(0 /H) = 8 .58), whic h has a sim ilar metallicity to local 
B-type stars dPeimbert 1987: Gies & Lambe rt 1992; Kilian 199i; 
iMathisllioOol) . For the LMC a nd the SMC, we adopt the oxygen 
abundance from lDufouJ ( Il984l) (12 + log(0/H) = 8.43 and 8.02, 
respectively). We confirm that metal-poor galaxies tend to have 
high dust temperatures. 

We should keep in mind that the story of star formation prop- 
erties in terms of metal licity is not so simple as we discussed above. 
lEngelbracht et al.l ( 1200 8) show that the dust temperature shows the 
peak around 12 + log(0 /H) ~ 8. Due to our small sample in such 
a low metallicity range, we cannot confirm it, but their results imply 
that the size and density of star-forming region is not a monotonic 
function of metallicity. There are also pieces of evidence that there 
are two categories of BCDs: "active" and "passive". The former 
hosts compact and intense star formati on activities, while the latter 
has d i ffuse star formation p roperties dHunt. Giovanardi." & Helod 
I2OO2I : IHirashita & HuriJ [20041) . Even with a similar metallicity, 
such a dichotomy exists, indicating that the star formation prop- 
erties are not determined simply by the metallicity. Active BCDs 
tend to be luminous in FIR because of large dust optical depth 
dHirashita & Huntll2004h . Thus, it would be possible that the cur- 
rent sample is biased to the "active " class. The presence of BCDs 
hosting cold dust in Virgo Cluster dPopescu et S]|2002l) also im- 
plies that such a bias could be present. This kind of bias should be 
examined in the future with more sensitive facilities. 

Finally, it is worth questioning which of dust-to-gas ratio and 
metallicity is more fundamental in regulating dust temperature. 
Larger absolute values of correlation coefficient found in Fig. |4| 
than those in Fig.|6|imply that dust-to-gas ratio rather than metallic- 
ity is more strongly connected with dust temperature. Thus, we here 
suggest that dust-to-gas ratio is more fundamental than metallicity 
in regulating the dust temperature. This should be further examined 
with a larger sample. 



5.4 Importance of submillimetre data 

It is expected that the radiative transfer effect, i.e. the shielding ef- 
fect of dust can be seen more clearly in submillimetre (submm) than 
in FIR, since the submm emissio n can trace dust with lower tem- 
perature (e.g. lGalliano et alj2003h . Although the submm data of the 
current BCD sample are still lacking except for those of IIZw40, 
future more sensitive submm facilities such as Herschel and ALMA 
will increase the data. 

We adopt A = 850 pm. as a representative submm wavelength, 
and examine the relation between (850/100)ci and (60/100)ci; 
that is, we adopt A = 850 pm instead of A = 140 pm, which 
is adopted in the rest of this paper. In Fig.|7] we present the colour- 
colour relation with Ay ~ (i.e. without the radiative transfer ef- 
fect; Section |231 . We also show the nea rby galaxy sample observed 
at A = 850 pm bv lDunne et"al] d2000!) . who also list the IRAS 60 
pm and 100 pm fluxes. We observe that the prediction traces the 
trend of the data, indicating that the submm flux is naturally ex- 
plained by the natural extension of the FIR SED. This means that 
a very cold component contributing only to the submm flux is not 
nece ssary for a m ajor part of the nearby galaxies. 

lGalliano"er al. (2003, 2005) argue that a very cold compo- 
nent should be introduced for metal-poor dwarf galaxies. In par- 
ticular, IIZw40, which is included also in the current sample, 
shows a clear excess of t he submm flux . In Fi g. |7l we also plot 
the data point taken from iGalliano et id] d2005l) for IIZw40 (see 
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Figure 6. Relation between FIR colour ({60/100)ci and (100/140)ci for Panels a and b, respectively) and oxygen abundance. The filled squares with vertical 
error bars represent the same BCD sample as in Fig.|4] The data points for the Milky Way, the LMC and the SMC are also presented by the open diamonds. 
The uncertainty in the oxygen abundance measurement is shown in Table|2] The correlation coefficients for the BCD data are r = —0.45 and r = 0.53 for 
Panels a and b, respectively. 



also lHunt, Bianchi. & Maiolino"2005^. Indeed, we obsei^e a clear 
deviation from the model prediction, which can be interpreted as 
a significant contribution from a very cold component. They also 
propose that such very cold dust should be located in an environ- 
ment where the stellar radiation is strongly shielded. In order to 
examine the effect of shielding (i.e. radiative transfer), we show 
the results with Ay > m Fig. [8] We observe that (850/100)ci 
increases even for Av > 2 while (60/100)ci is almost constant for 
such large Ay - This confirms that 850 /im flux is also sensitive to 
the shielded cold dust grains. Fig.[8]also demonstrates that the data 
point of II Zw 40 is explained with Ay ~ 5 with x ~ 300. Thus 
we conclude the existence of shielded cold dust in II Zw 40 (and 
in other galaxies wh ose data points deviate rightwards in Fig. [TJ. 
iGalliano et J] j2005[) also show that the UV ISRF of II Zw 40 is 
higher than that of the Milky Way by two orders of magnitude. 

Although our interpretation of the submm excess as a contri- 
bution from the shielded cold dust com ponent is in agreement with 
the analysis bv IGalliano et al. I j200 there is quantitati ve differ- 
ence between our results and theirs. Icralliano et al] ((2005) adopted 
Qabs oc A^^ at A > 100 /im. If the absorption efficiency is nor- 
malized to the value at A = 100 ^m, the absorption efficiency 
adopted in this paper (equation[3]l is 2 times larger than that used in 
lOalliano et alj (2005) at A = 850 iim. Therefore, w e require less 
very cold dust component than lOalliano et all ( 120051 ) by a factor of 
2. 



5.5 Common dust optical properties? 

In this paper, we have adopted a common FIR dust optical proper- 
ties; that is, we have applied the same absorption efficiency Qabs 
and grain size distribution as adopted in our previous model for 
the Milky Way and the Magellanic Clouds (HHS07). The same 
dust o ptical properties a re also consistent with the FIR colours of 
BCDs. lHibi et al.1 j2006h and HHS07 also argue that those dust op- 
tical properties also provide a good statistical fit to the FIR colours 
of nearby galaxies. Therefore, we suggest that the absorption ef- 
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Figure 7. (60/100)ci-(850/100)ci relations for silicate and graphite (as- 
terisks and open diamonds, respectively). The different points con'espond 
to X = 3, 10, 30, 100, and 300 from high to low (850/100)ci. The op- 
tically thin extreme for the ISRF (Ay = 0) is ad opted. Th e filled squares 
indicate the observational data of nearby galaxies in lDunne et al. (2000) and 
the op en triangle shows the data point of II Zw 40 taken from lGalliano et al] 
i2005l) . 



ficiency in the form of equation ([3]l, which iReach et"aD ( Il995h 
adopted to fit the dust emission spectra in the Milky Way, is gen- 
erally applicable in galactic environments. As shown above, the 
submm emission from nearby galaxies can also be explained by 
the same dust emissivity. 

It is als o important to stress that the classical dust emissiv- 
ity model bv lDraine&Leel ( ll984l) cannot explain the FIR colour- 
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Figure 8. (60/100)ci-(850/100)ci relations for (a) silicate and (b) graphite (asterisks and open diamonds, respectively). The results with x = 30 and 300 
are shown. For each value of x, Ay = 0, 0.5, 1, 2, 5, and 10 correspond to the data points from low to high (850/100 )ci (the point for v = 30 and Ay = 10 
is out of the range). Note that Ay = is the optically thin extreme for the ISRF. The nearby galaxy data taken from lDunne et alj i2OO0l) are presented by the 
filled squares. 



colour relation of the nearby galaxies teibi et al."200^; HHS07). 
As indicated by the emissivity assumed by Reach et al. ( 1995), it is 
better to adopt /3 ~ 1 for 100 /xm < A < Ai ~ 200 ^lm rather 
than to assume /3 ~ 2 in the entire FIR wavelength range longer 
than 100 /jm. Some amorphous materials indeed show /? < 2 (e.g. 
lAgladze et al.l ll996'). Amorphous grains are expected to form by 
the irradiation of cosmic rays ( Jiiger etal. 2003). 

Theoretically the break of the dust emissivity at Ai might be 
associated with the energy spl itting of the groun d state because of 
irregular amorphous structure. iMenv et aU l 20071) show that the ab- 
sorption coefficient of amorphous dust has a break at A ~ Am, 
which corresponds to the cut-off energy Tiujm of the energy split- 
ting due to the amorphous structure (tJm = 27r e/Am). Therefore, a 
possible interpretation is that Am ~ Ai. Since iMenv et alj j2007l) 
suggest Am ~ 700 /xm based on experimental data, the cut-off en- 
ergy corresponding to Am ~ 200 fim may be too high. Neverthe- 
less it is still worth investigating amorphous materials with higher 
cut-off energy as a candidate of cosmic dust. 



6 CONCLUSION 

We have investigated the properties of FIR emission of a sample 
of BCDs observed by AKARI, especially focusing on FIR colours 
and dust temperature. We have utilized the data at A = 65 fim, 
90 /im, and 140 /im, and have examined the relation between 
60 /xm-100 fim colour, (60/100)ci, and 140 /im-lOO nm colour 
(140/100)ci. Then, we have found that the FIR colours of the 
BCDs are located at a natural high-temperature extension of the 
DIRBE data of the Milky Way, the LMC and the SMC on the 
colour-colour diagram. We have explained the FIR colours also 
theoretically by assuming the same absorption efficiency, which 
may be appropriate for amorphous dust grains, and the same grain 
size distribution as the Milky Way dust. We have also shown that 
it is not easy to distinguish between a large dust optical depth and 



a low dust temperature only with FIR colours although addition of 
submillimetre data relax this degeneracy. 

In order to examine if the dust optical depth plays an impor- 
tant role in determining the dust temperature, we have investigated 
the correlation between FIR colour (dust temperature) and dust-to- 
gas ratio. We have found that the dust temperature tends to become 
high as the dust-to-gas ratio decreases as would be expected from 
the shielding effect of stellar radiation by dust. However, we fail 
to explain this trend quantitatively only by the effect of dust opti- 
cal depth. Thus, we conclude that there is a relation between dust- 
to-gas ratio and interstellar radiation field (ISRF) intensity. This 
relation is consistent with a "constant dust optical depth" of star- 
forming regions in BCDs, implying that dust extinction plays an 
important role in determining the condition for a burst of star for- 
mation. 

The correlation between dust-to-gas ratio and dust tempera- 
ture is e quivalent to that be t ween metallicity and dust temperature 
found in lEngelbracht et al. 1 1 I2OO8I) , since there is a correlation be- 
tween dust-to-gas ratio and metallicity. By comparing the correla- 
tion strengths, we propose that dust-to-gas ratio is more fundamen- 
tal than metallicity in regulating dust temperature, although we will 
have to confirm this with a larger sample. 
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